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Abstract 

In this talk we examine how one-loop soft and collinear splitting functions occur in the 
calculation of next-to-next-to-leading order (NNLO) corrections to production rates, and 
we present the one-loop gluon soft and splitting functions, computed to all orders in the 
dimensional regularization parameter e. Wc apply the one-loop gluon soft function to the 
calculation of the next-to-leading logarithmic corrections to the Lipatov vertex to all orders 
in e. 
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The single most important parameter of perturbative QCD is the strong coupling constant, 
a s , which has been determined in several ways [jjj. Some of the most promising ones are due to 
hadron production in e + e~ collisions; e.g., the hadronic branching ratio of the Z° or global event 
shape variables in e + e~ — ► 3 jets. The hadronic branching ratio Rz is known in perturbative 
QCD to three loops; however, the usefulness of this observable in the determination of a s is 
limited by the sensitivity of Rz to other Standard Model parameters [|| (for an overview, see 
ref. ||). On the contrary, e + e~ — ► 3 jets, which is known only to next-to-leading order (NLO) |4|, 
|5|, does not suffer from the above limitations. Thus a next-to-next-to-leading order (NNLO) 
calculation of this process could yield a significant reduction of the theoretical uncertainty in 
the determination of a s . 

In order to understand the general features of a calculation at NNLO, we begin by outlining 
how a higher-order calculation of a scattering process is performed. At leading order (LO) in 
a s the cross section is obtained by squaring the tree amplitudes. If n particles are produced 
in the scattering, each of them will be resolved in the final state. Thus no singularities appear 
in the LO cross section. At LO the coupling a s is evaluated with one-loop running, so that 
there is an implicit dependence on an arbitrary renormalization scale In addition, if one or 
both of the scattering particles are strongly interacting, the cross section will factorize into the 
convolution of parton density functions (to be determined experimentally) and a hard partonic 
cross section, which is computed as an expansion in a s . This procedure introduces into both the 
parton densities and the partonic cross section a dependence on a second arbitrary parameter, 
the factorization scale fj,p ||. Typically, the dependence on fiR and fip is maximal at LO. 

The calculation of the cross section at next-to-leading order (NLO) in a s is less straight- 
forward. Two series of amplitudes are required in the squared matrix elements: a) tree and 
one-loop amplitudes for the production of n particles; b) tree amplitudes for the production of 
n + 1 particles. The one-loop amplitudes typically have virtual ultraviolet and infrared singu- 
larities, which may be regularized using dimensional regularization. This involves analytically 
continuing the loop momenta into D = 4 — 2e dimensions, so that the one-loop amplitude is now 
a function of e. If this is expanded in e, the ultraviolet singularities appear as single poles in e, 
which can be removed by renormalizing the amplitude. This introduces an explicit dependence 
on the renormalization scale 

At NLO the structure of the infrared singularities has been extensively studied. Virtual 
infrared singularities appear as double poles in e, when they are both soft and collinear, and 
single poles in e, when they are either soft or collinear. Real infrared singularities occur in the 
phase-space integral over the n + 1 final-state particles of the squared tree amplitudes, either 
when any gluon becomes soft or when any two massless particles become collinear, thus yielding 
single poles in e. If one of the two collinear particles is soft, a double pole in e arises. The 
singularities occur in a universal way, i.e. independent of the particular amplitude considered. 
Accordingly, soft singularities can be accounted for by universal tree soft functions M, ||, and 
collinear singularities by universal tree splitting functions ||. These have also been combined 
into a single function [flOf] . A detailed discussion of the infrared singularities at NLO for e + e~ — > 



jets may be found, for example, in ref. [11]. 

For processes with no strongly interacting scattering particles, all infrared divergences cancel 
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when real and virtual contributions are put together to form the NLO coefficient in the expansion 
of the cross section Typically, the dependence on hr is reduced at NLO. For processes 

with strongly-interacting scattering particles, all infrared divergences cancel except for those 
associated with initial-state collinear singularities, which manifest themselves as single poles in 
e; these singularities are factorized into the parton densities, thus reducing the dependence of 
the cross section on fip ||. 

In order to compute a cross section at NNLO, three series of amplitudes are required: a) 
tree, one-loop, and two-loop amplitudes for the production of n particles; b) tree and one- 
loop amplitudes for the production of n + 1 particles; c) tree amplitudes for the production 



of n + 2 particles. For the case of NNLO e + e — > 3 jets the five-parton final-state tree [13 



amplitudes, as well as the four-parton final-state one- loop amplitudes exist in both helicity [14 



and squared matrix-element forms 15]. However, as we discuss below, in order to be used 
in NNLO computations higher-order terms in e must be included. For the required two-loop 
three-parton final-state amplitudes no computations exist, as yet. For single- and double-jet 



production at hadron colliders the six-parton tree [16, 17 1 amplitudes, as well as the five-parton 



one- loop amplitudes yi|, IS , 20] exist in helicity matrix-element form, but no four-parton two- 
loop amplitude computations exist, as yet. Indeed, no two-loop amplitude computations exist 
for cases containing more than a single kinematic variable, except in the special cases of maximal 



supersymmetry [21]. 



In the calculation of a production rate at NNLO the structure of the infrared singularities 
is the following: 

i) In the squared tree amplitudes, any two of the n + 2 final-state particles can be unresolved. 
Accordingly the ensuing soft singularities, collinear singularities, and mixed collinear /soft 
singularities have been accounted for by double-soft functions @ , double-splitting functions 



and soft-splitting functions [22], respectively. 



it) In the interference term between a two- loop amplitude for the production of n particles and 
its tree- level counterpart, all the produced particles are resolved in the final state and no 
new singularities appear through the phase-space integration. Thus, the expansion of the 
two- loop amplitude in e, which starts with a 1/e 4 pole, can be truncated at 0(e°). The 
universal structure of the coefficients of the 1/e 4 , 1/e 3 and 1/e 2 poles has been determined 

Hi ) In the interference term between a one- loop amplitude for the production of n + 1 particles 
and its tree-level counterpart any one of the produced particles can be unresolved in the 
final state; hence, the phase-space integration gives at most an additional double pole in e. 
Therefore, the expansion in e of the interference term starts with a 1/e 4 pole, from mixed 
virtual/real infrared singularities, and in order to evaluate it to O(e ), the (n + l)-parton 
one- loop amplitude needs to be evaluated to 0(e 2 ). (A similar need to evaluate one- 
loop amplitudes to higher orders in e has been previously noted in NNLO deep inelastic 
scattering [^4| and in the next-to-leading-logarithmic (NLL) corrections to the BFKL 
equation pi.) 
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iv) In the square of the one- loop amplitude for the production of n particles, the expansion 
in e of the amplitude, which starts with a 1/e 2 pole, must be known to 0(e 2 ) in order to 
evaluate the squared amplitude to O(e ). 

Here we focus on the singularities in in), which require that the (n + l)-parton one-loop 
amplitudes be evaluated to 0{e 2 ). For the case of NNLO corrections to e + e~ — ► 3 jets and to 
single- and double-jet production at hadron colliders, this would be a rather formidable task given 



the already non-trivial analytic structure of the one-loop e + e — > 4 partons amplitudes [14, 15 1 



and of the one-loop five-parton amplitudes [18, [T^, [2(J, both presented through O(e ) only. 
However, a simplification can be made if one uses the fact that the additional double poles in e 
of the interference term arise from the infrared-divergent regions of the phase-space integration. 
This implies that the one-loop (n+ l)-parton final-state amplitude needs be calculated to 0(e 2 ) 
only in the regions where two partons become collinear or one parton becomes soft. Therefore, 
one can use this amplitude calculated to 0(e ) and then supplement it in the soft or collinear 
regions by appropriate 0(e 2 ) terms. In these regions the amplitude factorizes into sums of 
products of n-parton final-state amplitudes multiplied by soft or collinear splitting functions. 
It is these soft or collinear splitting functions and the one-loop n-parton final-state amplitudes 
that must be evaluated to 0(e 2 ). This is a much simpler task than evaluating the full one-loop 
(n + l)-parton final-state amplitudes beyond 0(e ). 

Below, we provide the one-loop gluon splitting and soft functions to all orders in e [p9f|p|. 
A complete listing of the one-loop splitting and soft functions, including fermions, is given 
elsewhere Then we apply the framework outlined above to one of the effective vertices 

of the NLL corrections |31] to the BFKL equation ]32|] , namely to the one-loop amplitude for 



three-parton production in multi-Regge kinematics [33, 34 1, for which the produced partons 
are strongly ordered in rapidity. In NNLO and in NLL corrections to two-jet scattering, this 
amplitude appears in an interference term multiplied by its tree-level counterpart. Because of 
the rapidity ordering in the multi-Regge kinematics, the phase-space integration does not yield 
any collinear singularities; however, the gluon which is intermediate in rapidity can become soft. 
Accordingly the one-loop amplitude must be determined to 0(e ) plus the contribution with the 
soft intermediate gluon evaluated to 0(e) |25], [34]| . To determine the soft gluon contribution we 
use our all orders in e determination of the soft functions together with previous all orders in e 
determinations of the four-gluon amplitudes [35, 26, |36|]. 



We first briefly review properties of n-gluon scattering amplitudes. The tree-level color 
decomposition is (see e.g. ref . |37|] for details and normalizations) 



M* ree (l,2,...n) = g^-V £ Tr (TMDTMa) . . -T a 'W) m^ ee {u{l),a{2), . . . ,a(n)) , (1) 

a<£S n /Z 

n 

where S n /Z n is the set of all permutations, but with cyclic rotations removed. We have sup- 
pressed the dependence on the particle polarizations £j and momenta ki, but label each leg with 
the index i. The T ai are fundamental representation matrices for the Yang-Mills gauge group 
SU(N C ), normalized so that Tr (T a T b ) = 5 ab . The behavior of color-ordered tree amplitudes as 

* The one-loop splitting functions through O(e ) can be found in pj| , po| , and the one- loop soft functions 
through O(e ) may be extracted from the known four- [E7| and five-parton [FiL Esl E^] one-loop amplitudes. 
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the momenta of two color adjacent legs becomes collinear, is [37] 



m tree ^ ^ Split^ 6 (a Xa , b x * ) m*^ ( . . . if A . . . ) , (2) 

A=± 

where A represents the helicity, m^ ree are color-decomposed tree sub-amplitudes with a fixed 
ordering of legs and a and b are consecutive in the ordering, with k a = zK and fej = (1 — . 
For the case of only gluons, the tree splitting functions splitting into a positive helicity gluon 
(with the convention that all particles are outgoing) is 

Split t + ree (a + ,6+) = 0, Split^a" 6") ~' 



Vz(l-z) [ab] 

Spli1$-(a-&+) = SpHtf°(a+,6-)= (1 " f (3) 

V^(l - z) (a b) y/z(l - z) (a b) 



where the remaining ones may be obtained by parity. The spinor inner products pq , 17, 37| 
are (ij) = and [ij] = where are massless Weyl spinors of momentum fcj, 

labeled with the sign of the helicity. They are antisymmetric, with norm | (ij) | = | [ij] \ = y / s77, 
where Sjj = 2ki ■ kj. 

The behavior of color-ordered tree amplitudes in the soft limit is very similar to the above. 
As the momentum k of an external leg becomes soft the color-ordered amplitudes become 

r4 ree (..., a, fc± &, ...)|fc-o = Soft tree (a, fc± &) m^U-, a, 6, ...) , (4) 

with the tree-level soft functions 

The factorization of the collinear (^) and of the soft (Q) limits are similar. However, due to the 
locality of the collinear emission, the factorization property (0) extends to the full amplitude (||). 
Conversely, because of the non-locality of the soft emission and of the self-interactive nature of 
the gluon interaction, the factorization (||) is true only at the color-ordered amplitude level. 

The color decomposition of one-loop multi-gluon amplitudes with adjoint states circulating 
in the loop is p9[ 



Ln/2J+1 

M^(l,2,... n )= 5 " Yl E Gr n;i (a)m^ oop Kl),...,cT(n)), (6) 

where [x\ denotes the greatest integer less than or equal to x, Gr n; i(l) = N c Tr(T ai ■ ■ -T an ), 
Gr n;i (l) = Tr(T ai • • • T^'- 1 ) Tr(T a J • • -T a «) for j > 1, and S n; j is the subset of permutations 
S n that leaves the trace structure Gr n -j invariant, and where m*.' oop are color-decomposed one- 
loop sub-amplitudes. It turns out that at one-loop the m n; j>i can be expressed in terms of 
m n i°° P H0l> so we nee d only discuss this case here. The amplitudes with fundamental fermions 
in the loop contain only the m*.| oop color structures and are scaled by a relative factor of 1/N C . 
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The behavior of color-ordered one-loop amplitudes as the momenta of two color adjacent 
legs becomes collinear, is |2Cj 

m i-ioo P jsplit^V", b x >) m^°°P(. . . K x . . .) + Split ^ oop (a\ b Xb ) m^ x {. . . K x . . .)} . 

A=± ' 

(7) 

The one-loop splitting functions are, 

Split^ loop (a-,6-) 



(G f + G n )S P lit t + ree (a^,6- 
Split^ loop (a ± ,^) = G n SplitX ee (a ± ,b^), 



(8) 



Split^ loop (a + ,& + ) = -G f 



[a b] 



y/z{l-z) (ab) 2 ' 

The function G$ arises from the 'factorizing' contributions and the function G n arises from the 
'non-factorizing' ones described in ref. 0] and are given through O(e ) by H, |T| 



G n 



-If— 
e 2 Wl 



(9) 



7T 



) e + 21n(z)ln(l-z)-^- 

/ 



j(l - z){-s ab ) 
with Nf the number of quark flavors and 

1 r(i + e )r 2 (i-e) 



0(e) 



CT 



(10) 



(An) 2 -" T(l - 2e) 

As at tree-level, the remaining splitting functions can be obtained by parity. The explicit values 
were obtained by taking the limit of five-point amplitudes; the universality of these functions 
for an arbitrary number of legs was proven in ref. 41]. 

The functions (||) have been extended to all orders in e in ref. [^] 

(£)■,<!-,), (11) 



Gf 



2c r 



(3 - 2e)(2 - 2e)(l - 2e) 



G' 



cr 



1 



-Sab . 



l-Z 



7T€ 



+ 2 ]T e fc Li fc - 



where the polylogarithms are defined as |4 
Lii(«) =-ln(l-z) 

Lifc(^) = / — Li fe _i(i) 
Jo t 

and the regularization scheme parameter is, 

Sr = 



(fc = 2,3,...) 



£ 

n=l 



(12) 



1 HV or CDR scheme, 
FDH or DR scheme, 



(13) 



where CDR denotes the conventional dimensional regularization scheme, HV the 't Hooft- 
Veltman scheme, DR the dimensional reduction scheme, and FDH the 'four-dimensional helicity 
scheme. (For further discussions on scheme choices see refs. p7], S3].) 
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The behavior of one-loop amplitudes in the soft limit, as the momentum k of an external leg 
becomes soft, is given by 



m 



^ oop (...,a,fc ± ,&, ...)|fc— o 

l-loop 



(14) 



Soft tree (a, fc±, 6) m™i (-, a, &> •••) + Soft 1 - loo P(a, fc* 6) m^..., a, 6, ...) , 



where the one- loop gluon soft function may be extracted through 0(e°) from four- (27]] and 



five-parton |T^, Eq, 

Soft 1 - loop (a, fc* 6) = -Soft tree (a, fc ± , 6) c r 



one-loop amplitudes, and it is 



{-Sak){-Skb) 



i + y )+0(e). (15) 



Eq. (|l5|) does not depend on Nf or 5r. In ref. p9| we have extended it to all orders of e, with 
the result, 



Soft 



l-loop 



(0,^,6) 



-Soft* 



a, /c ± , b) cr 



7re 



1 / ^ 2 (-^) \ 

e 2 V {-s a k){-Skb)j sin(vre) ' 



(16) 



We now apply the results for the soft function ( |i~6| ) to the case of three-gluon production in 
multi-Regge kinematics. To do so, we also need the four-gluon one- loop amplitude through 0(e). 
In fact, this is known exactly to all orders of e. In the high-energy limit, s S> t, its dispersive 
part, which is all that contributes to the NLL BFKL kernel, is pjl 



Disp Ml- loop {A-,A' + ,B' + ,B-) = Ml iee (A-,A' + ,B'+,B-)g 2 c r I 



1 



2(1 - 2e) ip(l + e) - 2^(-e) + ^(1) + In ■ 



1 - ^ e ; 

H — - 4 

3 - 2e 



tj e(l - 2e) 

2(1 -e) 



(17) 



+ 



2e 



Nf 



where A, B and A', B' are respectively the incoming and outgoing gluons. The unrenormalized 
five-gluon one-loop amplitude in the multi-Regge kinematics, and in the soft limit for the inter- 
mediate gluon and to all orders in e, is obtained by using eq. (|14|), with the four-gluon one-loop 
amplitude ( |17| ) , and the dispersive part of the soft function ( |l6| ) , yielding 

DispMl- loop (A-,A>+, fc± B'+, B~)\ k _> = g 2 c r M 5 tree ( A~ , A'+ , fc± B'+, B")^ 



-5 + - U(l + e) - 2^(1 - e) + r/>(l) + In — 



(18) 



+ 



1-^e 
e(l - 2e) V 3 - 2e 



+ 



2(1 - e) 



e(l - 2e)(3 - 2e) 



Nf 



- N r 



(J, 



[1 + ei/)(l - e) - + e) 



which agrees through 0(e ) with the five-gluon one-loop amplitude, with strong rapidity ordering 
and in the soft limit for the intermediate gluon [18, [34|. Eq. (18) can than be matched to the 
full five-gluon one-loop amplitude, with strong rapidity ordering, computed through 0(e°). The 
result Q agrees with the NLL corrections to the Lipatov vertex computed in ref. |25[] in the 
CDR scheme, through 0(e). 
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In conclusion, in this talk we have examined how one- loop soft and collinear splitting func- 
tions occur in the calculation of NNLO corrections to production rates, and we have presented 
the one-loop gluon soft and splitting functions, computed to all orders in e. We have then ap- 
plied the one-loop gluon soft function to the calculation of the NLL corrections to the Lipatov 
vertex to all orders in e[E9|. A systematic discussion of the soft and collinear splitting functions, 



including the case of external fermions, is presented elsewhere [3C]. 
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